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The theory of stochastic transcription termination based on free-energy competitiontJ requires two 
or more reaction rates to be delicately balanced over a wide range of physical conditions. A large 
body of work on glasses and large molecules suggests that this should be impossible in such a large 
system in the absence of a new organizing principle of matter. We review the experimental literature 
of termination and find no evidence for such a principle but many troubling inconsistencies, most 
notably anomalous memory effects. These suggest that termination has a deterministic component 
and rnay conceivably be not stochastic at all. We find that a key experiment by Wilson and von 
Hippefl allegedly refuting deterministic termination was an incorrectly analyzed regulatory effect 
of Mg^+ binding. 

PACS numbers: 87.15.-v, 82.20.Pm, 61.43. Fs 
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I. INTRODUCTION 

The branching ratio of the termination process in gene 
transcription is balanced. In the case most thoroughly 
studied, p-independent termination in procaryotes, con- 
ventional gel experiments performed in vitro find a frac- 
tion P of elongating RNA polymerase reading through 
the termination sequence with | ln(l/P — 1)| < 4 essen- 
tially always, even though P is different for different ter- 
minators and can be made to exhibit order- 1 changes 
by perturbing the environment. This effect is astonish- 
ing from the standpoint of microscopic physics because a 
stochastic decision to read through or not requires a com- 
petition of transition rates - quantities of inverse time - 
that must be nearly equal for the branching to be bal- 
anced. RNA polymerase, however, is more the size of a 
glass simulation than a small molecule and thus possesses 
a broad spectrum of natural time scales spanning many 
decades. Without some physical reason for a particu- 
lar scale to be preferred, rate competition ought to have 
been severely unbalanced, meaning that one event occurs 
essentially always and the other never. Balanced branch- 
ing in termination has beien implicated in transcription 
regulation in a few cases,a but its broader significance, 
especially its robustness, is still a mystery. 

In this paper we examine the experimental facts rele- 
vant to the physical nature of termination with the goal 
of determining what, if any, principle selects the time 
scale for stochastic rate balance. Our conclusion is both 
surprising and unsettling. We find no evidence for such a 
principle, but glaring weaknesses in the case for stochas- 
ticity and a large body of unexplained experimental re- 
sults pointing to a termination decision that is partially 
deterministic. In light of the inaccessability of systems 
this large to ab-initio computation we conclude that tran- 
scription termination is a fundamentally unsolved prob- 



lem in mesoscopic physics and an ideal target for the 
emerging techniques of nanoscience. 



II. TERMINATION EFFICIENCY 

The simplest termination sequences are the p- 
independent terminators of procaryotes, which are capa- 
ble of causing polymerase to terminate in vitro in the ab- 
sence of the p protein factor. A representative sampling 
of these is reproduced in Table I. This diffops from lists 
that have appeared in the literature beforeaa by having 
been rechecked against the fuUy-sequenced genomeO and 
expunged of "theoretical" terminators identified only by 
computer search. They conform for the most part to the 
motif of a palindrome of typically 10 base pairs followed 
by a short poly-T stretch, although there is tremendous 
variety in the length and composition of the palindrome, 
variation in the length of the poly-T stretch, and occa- 
sional extension of the palindrome to include the poly- 
T stretch. This enormous variability contrasts with the 
simplicity of stop codons, which terminate protein syn- 
thesis by ribosomes and have no other function. 

p-independent terminators are characterized by "effi- 
ciencies" , i.e., the fraction of assayed transcripts that ter- 
minate. These rarely take on extreme values close to 1 or 
when measured in vitro. In cases where a measurement 
in vivo exists as well the latter is usually largerQ and is 
occasionally unmeasurably close to 1. Balanced termina- 
tion efficiency is commonly observed in vivo as well, how- 
ever. -Table II shows results from a particularly careful 
studyu in vitro in which termination probabilities in E. 
coli for wild-iypc terminators, mutant terminators, phage 
terminators Jj and terminators from S. Boydii were mea- 
sured under identical conditions. Despite the great vari- 
ety of these sequences the termination efficiency runs only 
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CGTTAATCCGCAAATAACGT AAAAACCCGC TTCG GCGGGTTTTT TTATGGGGGGA 
CAGTTTCACCTGATTTACGT AAAAACCCGC TTCG GCGGGTTTTT GCTTTTGGAGG 
CGTACCCCAGCCACATTA AAAAAGCTCGC TTCG GCGAGCTTTTT GCTTTTCTGCG 
ACACTAATCGAACCCGGCTCAAAG ACCCGC TGCG GCGGGT TTTTTTGTCTGTAAT 
AGTAATCTGAAGCAACGT AAAAAAACCCGCC CC GGCGGGTTTTTTT ATACCCGTA 
TCTCGCTTTGATG TAACAAAAAACCCCGCC CC GGCGGGGTTTTTTGTTA TCTGCT 
GAGTAAGGTTGCCATTT GCCCTCCGC TGCG GCGGGGGGC TTTTAACCGGGCAGGA 
CGATTGCCTTGTGAA GCCGGAGCGG GAGA CTGCTCCGGC TTTTTAGTATCTATTC 
CGTAAAGAAATCAGATACCC GCCCGC CTAATGA GCGGGC TTTTTTTTGAACAAAA 
GCGCAGTTAATCCCACA GCCGCCAG TTCCG CTGGCGGC ATTTTAACTTTCTTTAA 
AAATCAGGCTGAT GGCTGGTGACT TTTT AGTCACCAGCC TTTTTGCGCTGTAAGG 
AGGAAACACAG AAAAAAGCCCGCAC CTGACA GTGCGGGCTTTTTT TTTCGACCAA 
AGCACGCAGTCAAAC AAAAAACCCGCGC CATT GCGCGGGTTTTTT TATGCCCGAA 
CCCGTTGATCACCCATTCCC AGCCCCTC AATC GAGGGGCT TTTTTTTGCCCAGGC 
ACACGATTCC AAAACCCCGCCGG CGCAAA CCGGGCGGGGTTTT TCGTTTAAGCAC 
GAAACGGAAAACAGCGCCTG AAAGCCTCC CAGT GGAGGCTTT TTTTGTATGCGCG 
CTTAACGAACTAAG ACCCCCG C ACC GAAA GGT C CGGGGGT TTTTTTTGACCTTAA 
CCGCCCCTGCCAGAAATC ATCCTTA GCGAAACG TAAGGAT TTTTTTTATCTGAAA 
CATCAA ATAAAACAAAAGGC T CAGTC GGAA GACTG G GCCTTTTGTTTTAT CTGTT 
TCCGCCACTTATTAAGAAG CCTCGAG TTAACG CTCGAGG TTTTTTTTCGTCTGTA 
GCATCGCCAATGTAAATC CGGCCCGCC TAT GGCGGGCCG TTTTGTATGGAAACCA 
TGAATATTTTAGC CGCCCCAGTCA GTAA TGACTGGGGCG TTTTTTATTGGGCGAA 
ATTCAGTAAGCAGAA AGTCAAAAGCCTCCG AC CGGAGGCTTTTGACT ATTACTCA 
AGAAACAGCAAACAATCC AAAACGCCGC GTTCA GCGGCGTTTT TTCTGCTTTTCT 
CTGGCATAAGCCAGTTG AAAGAGGGAG CTAGT CTCCCTCTTT TCGTTTCAACGCC 
GCATCGCCAATGTAAATC CGGCCCGCC TAT GGCGGGCCG TTTTGTATGGAAACCA 
TGCGAAGACGAACAAT AAGGCC T CCC AAATC GGG G GGCCTT TTTTATTGATAACA 
ACGCATGAG AAAGCCCCCGGAAG ATCAC CTTCCGGGGGCTTT TTTATTGCGCGGT 
CATCAA ATAAAACGAAAGGC T CAGTC GAAA GACTG G GCCTTTCGTTTTAT CTGTT 
GGCATCAAATTAAGCAG AAGGCCATCC TGAC GGATGGCCTT TTTGCGTTTCTACA 
AAT TAATGTGAG TTAG CTCAC TC ATTA GGCACCCCAGGCTTTACACTTTATGCTT 
CTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTG 
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TABLE L p-independent terminators in E. coli taken primarily from Brendel et aljj These are oriented in the reading 
direction and are aligned at the poly-T stretch. The palindrome is underlined. Thef+)eginning and end of the selected sequences 
have no absolute meaning but simply follow the convention of d'Aubenton et alB The address identifies the location in the 
standard E. coli genomeH of the left-most nucleotide in the table. 



Sequence 


Name 


%T 


GGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAAT 


rrnB ti 

tonB t 

trp a 

trp a 1419 

trp a L126 

trp a L153 

T7Te 
T3Te 

tR2 
RNAI 


84 ± 1 
19 ± 1 
71 ± 2 

2 ± 1 
65 ±4 

8 ± 4 
88 ± 2 
14 ± 2 
49 ± 4 
73 ±4 


TCAAAAGCCTCCGACCGGAGGCTTTTGACTATTA 


CCAGCCCGCCTAATGAGCGGGCTTTTTTTTGAAC 


CCAGCCCGCCTAATGAGCGGGCTTTTGCAAGGTT 


CCAGCCCGCCTAATAAGCGGGCTTTTTTTTGAAC 


CCAGCCCGCCTAATAAGCGGACTTTTTTTTGAAC 


CTGGCTCACCTTCGGGTGGGCCTTTCTGCGTTTA 


GGCTCACCTTCACGGGTGAGCCTTTCTTCGTTCX 


GGCCTGCTGGTAATCGCAGGCCTTTTTATTTGGG 


AAACCACCGTTGGTTAGCGGTGGTTTTTTGTTTG 





TABLE II. Termination efficiencies measured m vitrou 
The first 3 terminators are native to E. c^li. These are fol- 
lowed by 3 mutants, 3 phage terminators^cl and one from 5*. 
Boydii. Far-right underlined sequences are termination zones. 



Sequence 



GTTAATAAC AGGCCTGC TGGTAATC GCAGGCCT TTTTATT 
GTTAATAAC AGGGGACG TGGTAATC CGTCCCC TTTTTATT 
TAATAAC AGGCCTGGC TGGTAATC GCCAGGCCT TTTTATT 
CCGGGTTAATAAC AGGCCTGC TTCG GCAGGCCT TTTTATT 
CGGGTTATTAACA GGCCTC TGGTAATC GAGGC TTTTTATT 
ATAACA GGGGACG TGGTAATC GCCAGCAGGCC TTTTTATT 
GTTAATAAAAGGCCTGCTGGTAATCGCAGGCCTTTTTATT 



GGTTCTTCTCGGCCTGCTGGTAATCGCAGGCCTTTTTATT 



Name 


%T 


tR2 


40 


tR2-6 


56 


tR2-ll 


54 


tR2-12 


69 


tR2-13 


11 


tR2-14 


20 


tR2-16 


36 


tR2-17 


67 



TABLE III. Termination efiieiencies for modified versions 
of the phage A terminator tR2.lij 



from 2% to 88%. Many other researchers report 
similar values for terminators in E. coli and othei; 
bacteria,EI including artificially altered terminators .1121 



Sequence 


Name 


rpo+ 


rpo203 


GCAACCGCTGGGGAATTCCCCAGTTTTCA 


trpC 301 
trpC 302 

trp t 
trp a 1419 
trp a 135 



8 

25 
3 
65 


20 
35 
45 
35 
80 


AACCGCTGGCCGGGATCGGCCAGTTTTCA 


CAGCCGCCAGTTCCGCTGGCGGCTTTTAA 


ACCAGCCCGCCTAATGAGCGGGCTTTTGC 


CAGCCCGCCTAATGAGCGGGCTGTTTTTT 





TABLE IV. Termination efficiences for wild- type p£'. coh 
polymerase (rpo+) and mutant polymerase (rpo203).t3 trp t 
is native to the genome. The rest are either mutants or syn- 
thetic. 

The results in Tables III and IV show balanced ter- 
mination for modified versions of the. phage terminator 
ti?lEll and for mutant polymeraseJl3 This also makes 
order- 1 changes to the efficiencies themselves.|_. Simi- 
lar effects were reported by other researcherdjO with 
different mutant polymerases. Modifications up to 20 
base pairs upstream and downstream of the terminator 
cause largfi-phanges to the efficiency without causing it to 
unbalance.13 Thus balanced termination efficiency is the 
norm rather than the exception. 



III. LARGE MOLECULES AND GLASSES 

Lacge systems are qualitatively different from small 
ones.til The specific heat of all non-crystalline matter in 
macroscopic quantities - including biokigical matter - is 
proportional to T at low temperatures.113 This behavior is 
fundamentally incompatible with the linear vibration of 
the atoms around sites, and is caused by collective quan- 
tum tunneling of atoms between energetically equivalent 
"frustrated" configurations Jl3 It contrasts sharply with 
the T^ behavior of crystals with small unit cells. Glasses 
also exhibit stretched-exponential time dependence in re- 
sponse to perturbations, i.e., of the form exp{~At^) with 
/3 < 1, indicating a broad spectrum of decay rates rather 
than just one. They also exhibit memory effects, such as 
"remanence" in spin glassesll3 or the well-known failure of 
ordinary silica to crystallize without annealing. This be- 
havior is universal and robust. All non-crystalline macro- 
scopic matter exhibits hysteresis, metastability, a broad 
spectrum of relaxation times, and memory. 

How large a system must be before it can exhibit 
such behavior is not known, as the relevant experiments 
are difficult to perform except on macroscopic samples, 
but there are many indications that even medium-sized 
proteins have glass-like properties. Crystals of myo- 
globin, a protein with a molecular weight of only 12,000, 
have linear specific heats at low temperaturesESl and 
exhibit stretched-exponential response to photodissoci- 
ation pulses .Ej Denatured proteins refold on a variety, 
of time scales ranging from nanoseconds to seconds,c2l 
and amino racids sequences chosen at random will not 
fold at alLtll Permanent misfolding of proteins with 
molecular weights of only 30,000 has been implicated in 
prion diseases.E3 Many enzymes exhibit hysteresis in their 



catalytic rates.l2jE3 The activity of cholesterol oxidase 
of Brevibacterium sp., a protein with molecular weight 
53,000, was recently shown by fiuorescence correlation 
techniques to have a memory effect persisting about j-i, 
second under normal conditions at room temperature.c3 
Other notable examples include wheat germ hexokkiase 
(mol. wt. 50,00CE3) with a half-life of 2 minutes£3 rat 
liver glucokinase (mol. wt. 52,00CE3) at 1 minute,t3 aud. 
yeast hexokinase (mol. wt. 50,000) at 1-2 minutes.c2l 
Thus RNA polymerase complexes, which have a molec- 
ular weight of 379,000 and are comparable in size to the 
largest computer simulations of glasses ever performed, 
are good candidates for systems that exhibit glassy be- 
havior. 

Classiness in enzymes is not always easy to observe. 
The mnemonic effect in yeast hexokinase occurs when it 
is preincubated with MgATP and free Mg^+ and the re- 
action is started with glucose, or preincubated with glu- 
cose and free Mg^+ and started with MgATP, but not if 
the enzyme is preincubated with glucose and metal-free 
ATP and then started with Mg^+.Ej Mnemonic behav- 
ior can be desjtepyed by "desensitizing" the enzyme with 
contaminants. c3 Time scales can depend on enzyme, sub- 
strate, product, activator and effector ligandrGOHceHtra- 
tions as well as pH, buffers, and temperature.c3E3'E2l Be- 
fore hysteresis and memory effects were recognized, early 
investigators generally adjusted such reaction conditions 
until the "improper" behavior was eliminated. E^ 



IV. POLYMERASE STATES 

While the size of RNA polymerase makes it plausible 
to expect glassy behavior on purely theoretical grounds, 
several direct lines of evidence indicate that the enzyme 
exhibits a spectrum of multiconformational, mnemonic 
and hysteretic behavior: 

1. Polymerase has a catalytic mode distinct from 
RNA synthesis, as it can cleave the RNA transcript 
through hydrolysis (rather than pyrophosphoroL: 
ysis, the reverse reaction of RNA syntkesis) ,EJ 
with the cleavage reaction requiring Mg^+ ,E2l being 
template-depeiident,E3 changing the polymerase 
footprint sizej^ asd— stimulated either by CreA, 
and CreB proteinskiea or by high pH (8.5-10.0).Ea 
The last effect was discovered serendipitously, go- 
ing unobserved for decades because assay condi- 
tions were being optimized to maximize elongatii 
rates, which occur at lower pH values (7.8-8.2tZl) 

2. RNA polymerase mobilities in non-denaturing elec- 
trophoresis gels show significant and discontinu- 
ous variance while bearing nearly identical tran- 
scripts or identical length transcripts with different 
sequences. l£| These mobility variances are still ob- 
served if the RNA traeBcript is first removed by 
ribonuclease digestion.La 



3. RNA polymerase ternary complexes vary greatly 
in their stability and mode of binding to DNA 
(ionic or non-ionic) in a template-dependent man- 
ner. Some complexes are stable against very high 
salt concentrations ([K+] — 1 M), while others 
(specifically those proximal to an upstream palin- 
drome sequence) are salt-sensitive (completely dis- 
sociating in concentrations as low as 20 mM K+). 
However, the salt-sensitive complexes axe stabilized 
by millimolar concentrations of Mg^+ o 

4. The size of the RNA polymerase footprint on the 
DNA template measured by ribonuclease digestion 
is significantly altered even at adjacent template 
positions, suggesting that the enzyme assumes dif- 
ferent conformations during elongation.cil 

5. Guanosine tetraphosphate (ppGpp) inhibits the 
rate of elongation on natural DNA templates but 
not on synthetic dinucleotide polymer templates, 
and does not inhibit elongation by competing with 
NTP binding, but by enhancing pausing. It must 
therefore bind to polymerase and modify its behav- 
ior at an unrelated regulatory site in an allosteric 
manner, rather than interfering with the substrate 
binding site.cHi 

6. The stability of a stalled elongation complex de- 
pends on whether the polymerase arrives art-ithe 
stall site via synthesis or pyrophosphorolysis.c^ 

7. Termination efficiencies are affected by transcribed 
upstream sequences and wntranscribed-downstream 
sequences adjacent to the terminator .cf 

8. Stalling elongating polymerase complexes (via nu- 
cleotide starvation) and then restarting them by 
nucleotide addition perturiss pausing patterns 50- 
60 base pairs downstream.E3 

9. An elongating polymerase's Michaelis constants Ks 
for NTPs vary cftffir 500-fold for different DNAptem- 
plate positions Jl3 and for different templates ,lI1 al- 
though these effects are not observed for synthetic 
dinucleotide polymer templates.cZl 

10. The rate of misincorporation at a single site 
for which the correct NTP is absent is signifi- 
cantly diffetent before and after isolation of ternary 
complexes .ES 

11. Stalled polymerase gradually "arrests" (i.e., is in- 
capable of elongating when supplied with NTPs), 
with the apppttximate half-tim&ipr arrest estimated 
at 5 minutesE^ and 10 minutesE^I for different DNA 
templates. The polymerase can continup-|elongat- 
ing if reactivated by pyrophosphorolysis.c3 

12. Even after undergoing arrest, crosslinking exper- 
iments show that the internal structure of poly- 
merase gradually changes over the course of the 
next hour.EJ 



13. Observations of single elongating RNA polymerase 
molecules show that it has two elongation modes 
with different intrinsic transcijiotion rates and 
propensities to pause and arrest .E3 

The possibility of metastability - through shape mem- 
ory or the conditional attachment of factors - is di- 
rectly relevant to the rate-balance conundrum because 
it provides a simple alternative to balanced stochastic 
branching that requires no physical miracles. If, for 
example, the polymerase possessed a small number of 
metastable configurational states and terminated deter- 
ministically depending on which state it was in, then 
balanced branching would be a simple, automatic con- 
sequence of scrambling the state populations. 



V. THERMAL ACTIVATION 

The idea that polymerase memory is potentially rel- 
evant to expression regulation is_iiot new.cZI It is im- 
plicit in tho-syork of Goliger et alEll and Telesnitsky and 
ChamberlinE3 and even explicitly speculated by the latter 
in print. However, because of the experimentaLevidence 
supporting the stochastic model of termination^ and the 
widespread belief - unjustified, in our view - that pro- 
teins equilibrate rapidly, this suggestion generated little 
enthusiasm. A key experiment supporting the stochastic 
model by Wilson and von HippelEI is both historically im- 
portant and typical, so it is appropriate that we consider 
it carefully. 

Wilson and von Hippel promoted and stalled RNA 
polymerase 8 base pairs upstream of the tR2 termina- 
tor hairpin of phage A in vitro, thermally equilibrated at 
temperature T, and then launched it forward by adding 
NTP. The results are reproduced in Fig. la. Termina- 
tion occurred at sites 7, 8, and 9 base pairs downstream 
of the beginning of the poly-T stretch (cf. Table II) with 
probabilities P7 = Nj/N, Pg = Nf^/N and Pg = N^/N. 
The data were originally reported as a semilogarithmic 
plot of 1/P — 1 against temperature, where P7 = Nj/N, 
Pg = Ns/iN - Nr) and Pg = Ng/{N - Nr - Ns). They 
concluded that all three branching probabilities P were 
thermally activated and had distinctly different activa- 
tion energies. However, it is clear from Fig. la that this 
conclusion is false. The three probabilities P are essen- 
tially the same function and are well characterized by the 
sum P = P7 -t- Pg -I- Pg , also plotted in Fig. la. This is 
shown more explicitly in Fig. lb, where the ratios P7/P, 
Pg/P, and Pg/P are plotted against temperature. The 
flatness of these curves shows that the branching ratios 
among the three sites are essentially constant and in- 
dependent of temperature within the error bars of the 
experiment. Note that these fractions are also all of or- 
der 1. Thus the alleged spread in activation energies was 
an artifact of the plotting procedure. 

Let us now consider the temperature dependence. It 
may be seen from Fig. la that P saturates to 1 at 80 
°C, the temperature at which Wilson and von Hippel 
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FIG. 1. a) Temperature dependence of termination proba- 
bility P forriphage A terminator tR2 reported by Wilson and 
von Hippel.a +, □, and x denote the probabilities to termi- 
nate 7, 8, and 9 nucleotides downstream from the beginning 
of the poly-T stretch. The sum is shown as O. b) -f , D, and 
X above divided by O to make a branching fraction /. c) 
Comparison of ionization model Eq. (1) with O from a). The 
ionization energy has been fit to eo = G.TeV (16 kcals/mole) 
and the quantity n/M'''''^ adjusted to make the curves match 
at 30 °C. d) Prediction of Eqn. (1) for dependence oh Mg ^ 
concentration compared with data of Reynolds et alE 

report that the polymerase "will not elongate", i.e., has 
stopped working properly. This suggests that the effect 
has something to do with the overall mechanical integrity 
of the enzyme rather than the termination process alone. 



Guided by this observation we note that the activated be- 
havior identified by Wilson and von Hippel is actually the 
formula for conventional monomolecular chemical equi- 
librium. The probability for a particle of mass M with a 
binding energy of Eq to be ionized off the polymerase is 



P = 



1 



1 + Ze^o/fciE 



"^?/. 



{\th = 



' 27r?i^ 
MksT 



(1) 



where n is the concentration of this component and Z 
is the change to the internal partition function that re- 
sults from binding. If one makes the approximation that 
\th is a slowly-varying function of temperature and can 
thus be taken to be constant then this reduces to the foig 
mula with which Wilson and von Hippel fit their data.EI 
That it works may be seen in Fig. Ic, where we plot the 
total termination probability from experiment against 
Eq. (1) with Eq — Q.leV and Z adjusted to match ex- 
periment at T = 30 °C. Thus reinterpreting this effect as 
an ionization equilibrium, we may account for the high- 
temperature intercept and weak temperature dependence 
seen in Fig. lb in the following way: In addition to the 
ionization state the polymerase possesses an internal con- 
figurational memory with a number of states of order 10. 
These code for termination at sites 7, 8 or 9. In the equili- 
bration step, the polymerase molecules come to thermal 
equilibrium and a fraction P of them become ionized. 
All of these terminate at one of the three sites when 
launched. The rest read through. 

A candidate for the ionizable component is an Mg^+ 
ion. In their studies of the effects of icm concentrations 
on termination efficiency, Reynolds et al3 discovered that 
Mg^+ has the strange and unique effect of increasing ter- 
mination efficiency to 100% for all terminators studied 
when reduced below 1 mM. The Mg^+ concentration in 
the experiments shovm in Fig. Id was 10 mM.0 Extrap- 
olating at T = 30 °CE3 using Eq. (1) we obtain, with no 
adjustable parameters, the fit to the [MgCl2] dependence 
found by Reynolds et afl shown in Fig. Id. The quality 
of this fit suggests that Mg^+ has a special function in 
regulating transcription, and that the temperature de- 
pendence in Fig. la is simply a thermal binding relation 
for this ion. This is cacroborated by the recent structural 
studies of Zhang et alj^ who report that polymerase crys- 
taUized out of 10 mM solution of MgCl2 has a Mg^+ ion 
bound at what appears to be the catalytic site of the 
enzyme. 

There is evidence for more termination channels other 
than the ionization ofJvIg^+. In Fig. 2 we reproduce re- 
sults of Reynolds ct alQ showing that terminator efficien- 
cies tend to saturate at large Mg^+ concentration to val- 
ues other than zero. The saturation values are balanced, 
and there is an evident tendency of them to cluster. Both 
effects are consistent with the polymerase executing an 
instruction at the terminator to read through condition- 
ally, even when the ionizable component is bound, if 
its memory is appropriately set. There is obviously not 




15 20 25 30 

[MgCls] (mM) 
FIG. 2. Termination efficiency as function of [MgCb] for 10 
terminators, as reported by Reynolds et al.Q The terminators 
are, top to bottom at the right edge, RNA I, T7Te, rrnB Tl, 
trp a L126, trp a, tR2, T3Te,P14, tonB t, and trp a L1.53. 

enough data here to draw such a conclusion, however. 
We note that Reynolds et alQ also found order-1 effects 
on the termination efficiency from Cl^ and K+, although 
with the opposite sign. The function of these ions is not 
yet known. 



VI. ANTITERMINATION 

What experiments can detect internal memory? In 
general, one would look for cases in which polymerase 
acts differently under apparently identical conditions, 
suggesting an internal control mechanism of some kind. 
Such thinking motivates the following hypothetical ex- 
periment: one constructs a template with promoter P 
followed by two identical terminators and flanking DNA 
sequences in succession. If termination is stochastic, then 
the branching ratio at T2 will be the same as that at Ti. 
If termination is deterministic and hysteretic, then the 
branching ratios will be different, depending on details. 
A passive termination at Ti would result in no termi- 
nation at T2, since the polymerase that reads through 
has been "polarized", i.e., selected for the memory set- 
ting that codes for read-through. An active termination 
at Tl would reprogram the memory there and cause a 
termination probability at T2 different from that of Ti 
but not necessarily zero. Variations of this design, e.g., 
adding more terminators, combining different termina- 
tors, changing their order, etc., could, in principle, an- 
swer more sophisticated questions, such as whether and 
how polymerase is reprogrammed in active read-through 
and whether non-equilibrium effects are important. 

A few such experiments have already been performed 
on DNA templates containing antiterminators (sequences 
upstream of terminators that reduce termination efficien- 
cies) and are thus less general than one would like, but 
they strongly support the idea of polymerase memory. 
There is indirect evidence in the case of N-antitermina- 
tion of phage A, the case most studied, that the memory 
is a physical attachment of the transcribed mRNA to the 



Sequence 


T7Te 


trp a 


AATTGTGAGCGGATAACAATTTCACACAGGAAACAGGGAA 


61 
51 
73 
45 
71 
75 
88 
99 


99 
52 
99 
99 
99 
66 
75 
80 


AATTGTGAGCGGATAACAATTTCACACAGGAAACAGAA . . 


AATTGTGAGCGGATAACAATTTCACACAGGAA . . . 


AATTGTGAGCGGATAACAATTTCACGGAA. . . 


AATTGTGAGCGGATAACAATTTCAGGAA . . . 


AATTGTGAGCGGATAACAATTTCGGAA . . . 


AATTGTGAGCGGATAGGAA . . . 


No Antiterminator 



TABLE V. Sequences and corresponding termination 
probabilities at downstream T7Te and trp a for modified|4ac 
antiterminators reported by Telesnitsky and Chamberlin.tJ 



Sequence 


oop t 


rpoC t 


AAATCTGATAATTTTGCCAATGTTGTACGGAATTC 


37 
45 
31 
29 
25 
17 
15 
11 
19 
20 


22 
17 
19 
16 
18 
20 
22 
20 
21 
16 


AAATCTGATAATTTTGCCAATGTTGGGAATTC . . . 


AAATCTGATAATTTTGCCAATGTTGGAATTC . . . 


AAATCTGATAATTTTGCCAATGGAATTC . . . 


AAATCTGATAATTTTGCCGGAATTC . . . 


AAATCTGATAATTTGGAATTC . . . 


AAATCTGATAATTGGAATTC . . . 


AAATCTGATAATGGAATTC . . . 


AAATCTGATAAGGAATTC . . . 


AAATCGGAATTC . . . 



TABLE VI> Antiterminator sequences constructed by 
Goliger et alEj from a promoter from phage 82, together with 
the readthrough probabilities in vitro for downstream termi- 
nators oop t and rpoC t. Note that these terminators are not 
in series. The underlined sequence on the right is the EcoRI 
linker. 

polymerase to form a lQQp.E3 There is also evidence that 
it is not true generally.^ i-. 

In 1989 Telesnitsky and Chamberlina reported mem- 
ory effects associated with the lac antiterminator found 
just downstream of the Ptac promoter in E. coli. Their 
key result is reproduced in Table V. Insertion of lac 353 
nucleotides upstream of the terminator makes different 
order-1 modifications to the termination efficiences of 
T7Te phage and trp a. The antiterminator contains a 
palindrome, and the antitermination effect is sensitive to 
modifications of the downstream 15-base-pair sequence. 
3 copies of T7Te placed in tandem downstream of lac 
showed that the antitermination effect is partially re- 
membered through multiple terminators: the efficiencies 
were 44%, 60%, and 90%, but without the antiterminator 
they were 90%, >90%, and >90%. 

In anothcE— .experiment in vitro reported in 1989, 
Goliger et alEll found that the E. coli terminator rpoC 
t and phage terminators oop t and tg,2 were strongly 
antiterminated by a sequence they constructed ac- 
cidentally. Their key result is reproduced in Ta- 
ble VI. A phage 82 promoter was fused onto a 
sequence containing either rpoC t alone or oop t 



Sequence 



GAGCGCGGCGG GTTCA GGA TGAAC GGCAATGCTGCTCATTAGC 
GCGTGGTCAAGGATGACTGTCAATGGTGCACGATAAAAACCCA 



Name 



putL 
putR 



TABLE VII. Antitermination sa 
from the Hong Kong phage HK022.E 



jiuences putL and putR 



and rpoC t in tandem using the EcoRI linker sequence 
GGAATTC. This resulted in unexpected antitermination 
in vitro of both terminators, but of different sizes that 
depended sensitively on the insertion point. The read- 
through effects in the tandem experiments were unfortu- 
nately poorly documented. One can see from Table V 
that the phage terminator responded more strongly in 
this experiment than did rpoC t. However, the reverse 
was the case in another experiment in which the antiter- 
minator was a portion of the 6S RNA gene downstream 
of a phage A pR' promoter, and in which factor NusA 
was present. As a control, this latter experiment was 
rerun with the phage terminator ^§2, which terminated 
at greater than 98% in all cases, seemingly immune to 
ant it erminat ioH,. 



King et alE3 reported in 1996 that the putL and 
putR antitermination sequences of the Hong Kong phage 
HK022,E3 shown in Table VII, caused downstream 
readthrough of a triple terminator consisting of tR' from 
phage A followed by the strong E. coli ribosome operon 
terminators rrn B <i and rrn B <2- This effect was sensi- 
tive to the choice of promoter. When putL was inserted 
between the Ptac promoter and the triple terminator 284 
nucleotides downstream and studied in vivo the termina- 
tion probability was 50%. Substituting the phage A Pl 
promoter for Ptac under the same conditions resulted 
in complete readthrough (though with wide error bars). 
When this experiment was repeated in vitro the antiter- 
mination effect was found to be smaller and to persist 
through all three terminators. The read-through proba- 
bilities at tR' were 34% and 31% for promotion by P^ 
and Ptac, respectively, but 57% and 27% for rrnB ti and 
76% and 40% for rrnB ^2- This result is incompatible 
with statistical termination, for both the antitermina- 
tion effect itself and the changes resulting from switching 
promoters are order-1 effects that do not add. They also 
reported that reduced Mg^+ concentration destroys the 
antitermination effect. 



VII. CONCLUSION 

In summary we find that the theory of stochastic ter- 
mination, which requires natural selection to engineer a 
physical miracle of balanced rates, is flawed, but that 
there is ample evidence of a sophisticated and as-yet 
poorly understood regulatory system in RNA polymerase 
involving hysteresis, metastability, and long-term config- 
urational memory, all robust phenomena in inanimate 
matter. On this basis we predict that branching ratios of 



identical terminators in series will differ by order-1 
amounts very generally - specifically in the absence of 
looping. We propose that the confusion surrounding 
the existence of polymerase memory is symptomatic of 
the larger problem that measurement of physical activ- 
ity on the length and time scales appropriate to life has 
thus far been impossible, and that overcoming this prob- 
lem should be one of the high-priority goals of modern 
nanoscience. 
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